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Abstract. In ion traps, entangling gate operations can be realized by a bichromatic 
pair of laser beams that collectively interact with the ions. In this paper, a new method 
of modelling the laser-ion interaction is introduced that turns out to be superior to 
standard techniques for the description of gate operations on optical qubits. The 
treatment allows for a comparison of the performance of gates based on ® fz and 
<y<l> ® C0 interactions on optical transitions where the bichromatic laser field can be 
realized by an amplitude-modulated laser resonant with the qubit transition. Shaping 
the amplitude of the bichromatic laser pulse is shown to make the gates more robust 
against experimental imperfections. 



1. Introduction 

The processing of information based on the laws of quantum physics [T] has become a 
very active field of research during the last decade. For the experimental demonstration 
of fundamental key results of quantum information theory, ion-trap based systems 
have played a major role. The success of ion trap experiments can be attributed 
to the fact that encoding quantum information in either hyperfine ground or meta- 
stable excited atomic states provides well-defined quantum bits (qubits) with long 
coherence times. The use of lasers for manipulating the qubit state allows for precisely 
switchable interactions with low decoherence rates. The fundamental operations of 
qubit initialization, arbitrary single qubit manipulation and quantum state-detection 
have already been used in atomic clocks with single ions for many years. Contrary to 
other realizations [1] of quantum information processing, the most demanding operation 
in ion traps consists in the realization of an entangling gate operation. Because of the 
repulsive Coulomb force, the inter-ion distance is orders of magnitude bigger than the 
characteristic length scale of any state-dependent interaction between ions in ground or 
low-lying excited states. In all current experiments creating entangled ions [21 [3l HI [5] , 
gate operations rely on interactions that are mediated by the vibrational degrees of 
freedom of the ion string. These gate operations fall into two categories: 
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(i) Quantum gates induced by a laser beam that interacts with a single ion at a time 
as originally proposed in the seminal paper by I. Cirac and P. ZoUer [6] and later 
realized by the Innsbruck ion trapping group [2J. In these gates, a single ion is 
entangled with a vibrational mode ^ of the ion string and the entanglement is 
subsequently transferred from the vibrational mode to the internal state of a second 
ion. 

(ii) Quantum gates induced by a bichromatic laser that collectively interacts with two 
or more ions. Here, a vibrational mode becomes transiently entangled with the 
qubits before getting disentangled at the end of the gate operation, resulting in an 
effective interactions between the qubits capable of entangling them. Gates of this 
type were first proposed by G. Milburn P[9], A. S0rensen, K. M0lmer [IDIIII] and 
E. Solano [12], and subsequently realized by ion trapping groups in Boulder, Ann 
Arbor and Oxford [S, U [5] . 

Even though both classes of gates are applicable to hyperfine qubits as well as optical 
qubits (i.e. qubits encoded in hyperfine states or in states linked by a dipole-forbidden 
transition with an optical wavelength), current experiments with optical qubits have 
relied on the former and experiments with hyperfine qubits on the latter type of 
interaction. In any case, the main goal consists in demonstrating fast operations creating 
entanglement with high fidelity. 

The purpose of the present paper is to discuss bichromatic gate operations with 
a focus on implementations using an optical transition. It turns out that for optical 
transitions, gate operations are achievable by illuminating the ions with an amplitude- 
modulated laser beam that is resonant with the qubit transition. The paper is organized 
as follows: section [2] reviews different methods of realizing bichromatic quantum gates 
and discusses properties that are specific to their application to optical qubits. In 
section [3l an effective Hamiltonian for the laser-ion interaction will be derived by going 
into a reference frame rotating at non-uniform speed in order to eliminate non-resonant 
excitations of the qubit transition that do not couple to the vibrational mode. In this 
way, it will be shown that for a single ion qubit the interaction is well described by a 
Hamiltonian H = ih{'y{t)a^ — 'y{t)*a)a^ where the coupling strength 7 is proportional 
to the laser intensity in the limit of low intensities but starts to saturate at higher 
intensities and where = a ■ is a. component of the Pauli spin operator a coupling 
to a vibrational mode of the ion described by creation and annihilation operators a\ 
a. Furthermore, it will be shown that depends not only on the particular type of 
gate operation but also on the laser intensity and the relative phase between the two 
frequencies of the bichromatic field. Equations (|TOl) . (ITTj) . ( 1271) describing the action of 
the gates based on (S> cr^ and on o"^ a,^ (M0lmer-S0rensen) interactions are the key 
results of the paper. For the M0lmer-S0rensen gate, the result will be compared to the 
analysis presented in [11]. In addition, the performance of <7z®<Jz and a^^a^ gates will 
be compared. Section H] shows how to use pulse-shaping of the laser intensity as well as 
spin echo techniques to make the gates more robust against fluctuations of the control 
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parameters. 

2. Quantum gate operations based on bichromatic laser fields 

2.1. Driven quantum harmonic oscillator 

The Hamiltonian H = hva'^a + hrLi{a'^e^^^ — ae"*"^*) describes a liarmonic oscillator 
oscillating at frequency v and driven by a force with frequency uj and coupling 
strength Vt. Going into an interaction picture defined by Hq = hva^a yields the 
Hamiltonian H = hQi{a^e^^^ — ae~^^^), where 6 = u — u. Under the action of the 
driving force, an oscillator that is initially in a coherent state remains in a coherent 
state. For a force that is slightly detuned from resonance, the coherent state maps out 
a circle in phase space and returns to the initial state after a period r = 27i/6. This 
operation multiplies the oscillator state by a phase factor whose magnitude is given by 
the ratio of the strength of the force and the detuning as shown in p5]. 

In order to allow for variations of the driving field's strength, we generalize the 
Hamiltonian to H = ih{-y(t)a^ — 7*(t)a) and calculate its propagator U{t) by using 
the Baker-Campbell-HausdorfF relation D{a)D{f3) = D{a + (3) exp{ilm{aj3*)) for the 
displacement operator D{oi) = e°"^~"*". For the propagator, we find 

n 

U{t) = lim Y\exp{-lH(tk)At) = D{a{t)) exp{t^t)) (1) 

k=l 

where At = t/n, tk = kAt and 

a{t) = [ dt'^{t'), 
Jo 

<^{t) = lm [ dt'-i(t') [ dt"-i*(t"). 
Jo Jo 

In case of a driving force with constant amplitude, 7(t) = Qe'^^^, one obtains a(t) = 
I (f) (1 - e'^') and $ = (f)^ {5t - sm5t). After a time tn = 27rN/\S\, N = 1,2,..., 
the coherent state returns to its initial state in phase space with its phase changed by 
an amount $(rAr) = 2ttN (^) sign (5). By making this phase change depend on the 
internal states of a pair of ions, an entangling gate operation can be achieved. For 

H = th{^{t)a^ -^*{t)a)0 (2) 

where O is an operator acting on the qubit states, the propagator ([T]) is replaced by 

U^{t) = b{a{t)0) exp(i$(t)02). (3) 

Choosing the interaction time r such that «(t) = thus realizes a propagator that 
depends nonlinearly on O and does not alter the vibrational state. 
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2.2. Laser-ion interaction 

The interaction of a single ion qubit resonantly excited by a monochromatic laser field 
with frequency uol is usually described by performing a rotating- wave approximation 
with respect to the optical frequency to obtain the Hamiltonian 



Here, 6 = ool — ujq is the detuning from the qubit transition frequency ujq, v is the 
frequency of the ion's vibrational mode of interest, and cr+ = + i^y)l'^ with the 
Pauli matrices (Jx,y The strength of the laser- ion coupling is characterized by the Rabi 
frequency f2, and the strength of processes involving changes in the vibrational state 
is determined by the value of the Lamb-Dicke parameter rj. Equation (jll) represents a 
Hamiltonian in an interaction picture that is defined with respect to the Hamiltonian 
Hq = hua'^a-^- describing the ion qubit in the absence of any laser-ion interactions. 
If ^ 1, the Lamb-Dicke approximation e"^^^""^ ""'+ate''") ^ 1 -|-ir^(ae~*^* + a^e**^*) is used 
to simplify (jl]). The resulting three terms describe excitations on the carrier, the lower 
and the upper motional sideband, respectively. The generalization of the Hamiltonian to 
the case of two and more ions is straightforward. For the sake of simplicity, calculations 
in section [3] will be limited to the case of the laser coupling to the centre-of-mass mode 
along the axis of the ion string where all ions experience the same coupling strength. A 
detailed account of laser- ion interactions is given in [T3] . 

2. 3. Gz ® cFz gate 

A Hamiltonian as described by ([2]) was employed for the first time in an experiment 
[H] creating a Schrodinger cat state with a single ion using O = az, i.e. a coupling to 
the motional mode that depended on the internal energy eigenstate of the ion. Later, 
it was realized that the same type of coupling could be used to entangle a pair of ions 
by performing a conditional phase gate P, |3]. In the experimental realizations [31 [5], 
spin-dependent forces acting on a pair of hyperfine or Zeeman ground states have been 
realized by near resonant driving of Raman transitions between vibrational states (see 
Figured]). For this purpose, two non-copropagating laser beams with frequencies uji,,ujr 
form a moving standing wave with difference frequency ujb — uJt close to the frequency 
z/ of a vibrational mode. The ac-Stark shift of the qubit states | |),| t) results from a 
non- resonant coupling to another atomic state |e) that is made qubit state-dependent 
by properly chosen polarizations. Since the laser field exhibits a strong spatiotemporal 
modulation, the resulting potential gradients induce a force acting on the qubits that 
is state-dependent and that couples to the vibrational mode by displacing the qubit 
along a circle in phase space. When the ions couple to the centre-of-mass mode (stretch 
mode), the coupling to the mode can be made to disappear when both qubits are in the 
same quantum state by choosing an ion spacing that is an odd (even) integer multiple 
of half the wavelength of the moving standing wave. Disregarding unimportant global 
and single qubit phases, this coupling is then described by an operator O = Sz where 




) + h.c.. 



(4) 
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Figure 1. Raman coupling of motional states for a hyperfine or Zeeman qubit. (a) 
The states \ ] — Q) ,n — 1) are coupled via the excited state |e). A similar coupling 
not shown in the figure exists for the qubit state | [). The laser detunings Ai, A2 from 
the mediating states are large compared with trap frequency v and the qubit level 
spacing to avoid spontaneous emission from state |e). Constructive interference 
of paths 1 and 2 is achieved for counter-propagating laser beams, (b) The coupling 
is maximized for counter-propagating laser beams with frequencies Ub, oJr forming a 
moving standing wave along the axis of vibration of the motional mode of interest. 



Sz = (Tz + Oz i?> a. collective spin component of the qubits. 

The situation is different for qubits encoded in atomic states connected by a narrow 
optical transition. For coupling motional states | |, n = 0) ^ | |, n = 1), here, the 
other qubit state | f) serves to mediate the coupling. Similarly, a coupling between 
the states | 1,72 = 0) ^ | t,""- = 1) is mediated by the state | [). To achieve a strong 
coupling, a detuning Aj from the intermediate states can be chosen that is smaller than 
the transition frequency between vibrational states provided that the decay rate of the 
metastable state is small compared to v. For ujh-r = cuq ± z//2, the two interfering paths 
shown in Figure [2^ connecting levels 1 1, n = 0) ^ 1 1, n = 1) have equal strength. Since 
the detunings from the mediating states now have opposite signs, destructive interference 
is achieved for counter-propagating beams whereas the coupling is maximized for co- 
propagating beams. In the limit of small excitation (fi <^ z/), the coupling strength VLr q 
on the Raman transition between | |,n = 0) and | |,n = 1) is given by = 2riD? /v. 
The states | t,n = 0) and | t,n = 1) are coupled with equal strength but opposite sign. 
For stronger excitation, the carrier transition is non-resonantly excited which leads to 
a saturation of VLr^. As long as the intensities of the bichromatic beams are equal, 
there is no overall ac-Stark shift due to excitation of the carrier transition and the first 
motional sidebands because ac-Stark shifts caused by the two laser fields exactly cancel 
each other. 

In the case of two ions excited on the centre-of-mass mode, a driven quantum 
mechanical oscillator is realized with collective atomic oscillator O = Sz- In addition 
to the coupling of vibrational states, there is another small M0lmer-S0rensen coupling 
[10] that does not exist for the case of hyperfine or Zeeman qubits: collective spin flips 
between the states | [[,n) and | TT;''^) occur by processes involving a blue and a red 
photon that are mediated by the states | J,t,r2 ± 1) and | |J,,n ± 1) (see Figure [3]). A 
similar process involving either two blue or two red photons couples the states | Ti;^) 
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Figure 2. Raman coupling of motional states for an optical qubit. (a) The coupling 
for states \ ^ ,n = Q) ,\'\ ,n = \) is mediated by the states | i, n = 0), 1 1, n = 1) and vice 
versa on the narrow qubit transition. Since spontaneous scattering from the mediating 
state is small, the detuning can be made small compared to the trap frequency. For 
Ai = A2 = —1^/2, the coupling is maximized by choosing a copropagating beam 
geometry, (b) Optimum coupling is achieved for co-propagating laser beams with 
frequencies ojf,, Ur propagating along the axis of vibration of the motional mode of 
interest. As the bichromatic laser field could also be described by a monochromatic 
laser field that is amplitude-modulated with a frequency LOt — LVr close to the vibrational 
frequency, no spatially varying ac-Stark shifts are involved in the coupling. 




Figure 3. M0lmer-S0rensen gate. A bichromatic laser field with frequencies U!b,ujr 
satisfying 2wo — ujb + is tuned close to the upper and lower motional sideband of 
the qubit transition. The field couples the qubit states | H) | tt) via the four 
interfering paths shown in the figure. Similar processes couple the states | ti) | it) 
with the same strength provided that the Rabi frequencies of the light fields Ub, cUr sue 
equal. 



and I it, n). 

2.4- (^4,® (^4, gate 

In contrast to ® gates that do not change the internal states of the ions, the 
C(/> ® C(/> gate operations first investigated by A. S0rensen, K. M0lmer [10] and others 
[12] rehes on collective spin flips | ||) | It)) I it) ^ I it) by processes coupling 
to the lower and upper motional sidebands as illustrated in Figure [3l It can be 
shown that the Hamiltonian governing the action of the gate is described by setting 
O = coscpSx + sin (j)Sy [151 HB]- For a properly chosen coupling strength, the gate 
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operations maps the product state basis {ITT)?!!!))!!!)?!!!)} oiito a basis of entangled 
states. For hyperfine qubits, a detailed discussion of advantageous beam geometries is 
presented in |18]. In the case of optical qubits, it is again possible to choose a pair 
of co-propagating beams for performing the gate operation. The only difference to the 
(^z^^^z gate consists in the choice of laser frequencies ujb^r = ojo^p required for achieving 
a resonant coupling. Formally, the gate operation is equivalent to a cr^ (S> Uz interaction 
in a rotated basis. To stress this analogy, the M0lmer-S0rensen gate operation is often 
also called & a^® gate. 

The possibility of choosing co-propagating laser beams for performing either cr^^a^ 
or o"0 ® gates is attractive from an experimental point of view. The light field could 
be generated by passing a laser beam through an acousto-optical modulator driven by 
two radio-frequency fields and subsequently coupling the first-order diffracted beams 
into a single-mode optical fibre, thus reahzing a simple and stable setup 0. If the Rabi 
frequencies f^b, of the blue and the red detuned laser beam are equal, light shifts 
due to the non-resonant excitation of the carrier transition and the first-order sidebands 
are exactly cancelled. Light shifts arising from coupling to other Zeeman transitions or 
far-detuned dipole transitions could be cancelled by suitably balancing the ratio VLi,/VLt. 

3. Effective Hamiltonians for o"^ ® and ® gates 

We are interested in deriving an effective Hamiltonian that accurately describes the 
dynamics on a pair of optical qubits induced by a co-propagating bichromatic laser field 
with frequencies loi^^r = i^o =t where the detuning is either close to half the vibrational 
frequency or close to the vibrational frequency, i. e. b = {v—t)/2 oi 5 = v—t with t <^v. 
As optical qubits interacting with lasers typically have smaller Lamb-Dicke parameters 
than hyperfine qubits coupled by Raman transitions, transient non-resonant excitation 
of the carrier transition is expected to play an important role for < z/. Whereas 
usually non-resonant interactions are taken into account only qualitatively after having 
derived an effective Hamiltonian, in the following calculation they will be eliminated 
right at the beginning by going into a reference frame rotating at non-uniform speed. 
The Hamiltonian for the bichromatic laser field we are interested in, is given by 

The laser is assumed to interact collectively with m ions on the axial centre-of-mass 
mode. Here, S+ = Ym=i'^+'^ '^^ denote operators annihilating and creating 

phonons. It is also possible to interpret this interaction as being due to a single resonant 
laser beam that is amplitude-modulated with modulation frequency 6. The optical phase 

I In principle, cr^ (g) cr^ gate operations could also be achieved by phase-modulating a laser with a 
modulation frequency close to v by choosing a modulation index where the carrier strength vanishes. 
However, this approach has the strong disadvantage that even small changes of the modulation index 
from the desired value give rise to light resonant with the transition which is having desastrous effects 
on the gate performance. Similarly, modulating at frequency close to 2v and using the carrier and one 
of the first sidebands is problematic because of light shifts induced by the other sidebands. 
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of the laser field is denoted 0, and the phase C, accounts for a time difference between the 
start of the gate operation and the maximum of the amplitude modulation on the laser 
beam. Using the picture of an amplitude-modulated resonant beam, it is obvious that 
there are fast dynamical processes on the carrier transition with a periodicity given by 
r = 27r/(5 that excite the ions to the other state in the first half of the period and transfer 
it back to the original state in the second half. We are not really interested in exactly 
calculating the dynamical evolution of the quantum state on this fast time scale. Rather, 
we would like to know the time evolution at the instances r, 2r, 3r, ... ,t = 271/5. It is 
useful to rewrite ([5]) as 

H = hf{t){e~''^S+b{i7]e"'') + e"''S-i){-irie"'')) 

= hf(t){{Sx cos (p+Sy sin (t>){D^ + D^) + i{Sy cos — 5*0; sin 0)(-D+ — -D_)) 

= : + D^) + ^Sl^\D^ - D-)) , 

where we used the displacement operator D{a) = and the definitions D± = 

D{±ir]e'''')/2, Si'''^ = 5^ cos + sin 0, si'^^ = 5^ cos - 5^. sin and f{t) = 
2flcos{St + C). For 6 = 1^/2 or 6 = u as required by either Cz^cTz or (j<^®cr<^ interactions, 
the Hamiltonian is periodic in time, i. e. H{t + r) = H(t) with period r = 271/6. Note 
that we assume the two-photon couplings to be strictly resonant for the moment. In 
the next step, we will get rid of the fast non-resonant carrier oscillation by going into 
another interaction picture defined by Hq = hf{t)Si'^\ Writing 

H =hf{t)S'^f + Hi, with 

H, = hf{t){Si^\D^ + D_-1)+ tSlf\D^ - D_)) , 
we obtain the interaction Hamiltonian 

= hf{t)Si'^\D+ + D^-l) + hfit) {cos{2F{t))Sf - sin(2F(t))5,) i{D+-D^), 
where 

2Q 

F(t) = — (sin(5t + C)-sinC). (6) 

Now, we can approximate the time evolution over the course of an oscillation period by 
a Magnus expansion of the propagator [19j in order to obtain 



Uj{t) = exp I -1 dt'H,{t') -^1^ dt' 1^ dt"[Hj{t'), Hj{t")] + • • I ■ (7) 

From now on, the phase will be set to zero to simplify the notation. The formulas in 
the remainder of section [3] are easily generalized to the case of arbitrary by making 
the replacements Sx — ^ S^"* and Sy — > S\f'\ 

In the following two subsections, effective Hamiltonians for the cr^ ® gate and 
the ® o"0 gate will be derived starting from ([7]). We are interested in obtaining 
Hamiltonians of the form given in ([2]) that are valid in the regime f2 ^ z/. In addition, 
the calculation is going to yield correction terms for the case when ^ z/ no longer 
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strictly holds and additional terms that do not commute with the atomic operator O in 
([2]). Towards this aim, terms proportional to higher orders of the expansion parameter 
will be dropped. In the calculation, Bessel functions J„(x) will be evaluated at 
X = 4Q/6. These functions will be kept till the end of the calculation and expanded in 
Q/6 only for the final analysis. 

3.1. (Jz® cFz gO'te 

For the ® gate, we set 5 = v jl. We start by calculating the first term 
a'^^ff ~ \ jo appearing in the exponent of ([7]). Here, it is important to 

note that the integrand Hj{t) = YlkL-oo ^(k)^^^^^ ^ periodic function of time with 
period r = 2tt/6 so that all its Fourier components except the constant term iif(o) will 
average to zero when they are integrated over one period. The only non-zero Fourier 
components of the function f{t) = fi(e^('^*+f) +e-*('^*+^)) = J2n=-oo fne'""^^ are /+i = Qe'^ 
and /-I = i^e"*''. Since all non-zero Fourier components of (D^+D^—l) are even integer 
multiples of 6, the S'^^-term of Hj averages to zero. Therefore, we obtain 

= - f dtfit) {cos{2F{t))Sy - sin(2F(t))^,) i{D+ - D.) . 
Jo 

In the Lamb-Dicke limit, 

oo 



AnSt 



and the components (i+2 = — ?7a^ and (i_2 = —rja are the only relevant ones. Finally, we 
have 



^gi^sm(<5t+C) _^ ^tg-i^sin(<5t+C) 
o 



n=— OO 
oo 

AnSt 



n=— oo 



where J„ is a Bessel function, A = ^{Sy + iSz)e ''^ with 
AVt ^ 

V' = ^smC, (8) 

and an = {AJn{f) + At J„(-^))e*"^ = {A+ (-l)"At) J„(^)e^"f. In the following, the 
argument 4:Q/6 of the Bessel functions J„ will often be dropped to keep the notation 
simple. It is convenient to express A± A^ a.s 

Sy^^ := Sy COS 'ip + SzSin-ip = A + 

Sz,xi) '■= Sz cos^jj — Sy simp = —i{A — A^). 



Ion trap quantum gates with amplitude-modulated laser beams 



10 



Note that the hnear transformation ^ preserves the usual Lie algebra commutation 
relations for the operators S^, Sy^^, S^^^. The four terms /+i(i„2'3+i,/-i'i+2'3-i, and to a 
lesser degree /+i(i4.2a-3, /_i(i_20.+3, contribute to H^ff Evaluating 

/+id_2a+i = (fie^^)(-r/a)(A - A^)J,{^)e< = -r^nj^e^'^iaS,,^ 
as well as the other terms, we arrive at the effective Hamiltonian 

where r]Q{Ji + J3) ^ [27]^"^ /6){1 — 4i7^/(35^)). This Hamiltonian describes a spin- 
dependent force that starts to saturate when the Rabi frequency goes up. While the 
atomic operator O = Sz,ip coincides in the limit of weak excitation with the operator 5*2 
obtained from second-order perturbation theory, it depends on the phase C between the 
blue- and the red-detuned laser beams in the limit of strong excitation. For the periodic 
Hamiltonian Hi(t) = J2T=-oo^k€-^^^* , the second order contribution to the effective 
Hamiltonian H^ff given by 

1 °° 1 

m=l 

After evaluating the commutators [if(3), if(_3)], the effective Hamiltonian 

is obtained (the contribution of the commutator [H(2),H(-2)\ oc {rjVLY {VL / 5)^ is 
insignificant). If the detuning 5 = (z/ — e)/2 slightly deviates from half the oscillation 
frequency the Hamiltonian is given by 

= ihr^n{J, + J3)^.,^(ate^(^*-2C) - ae'^^^'-^^)) - Jo^^J,^ (10) 

is obtained. The second order term H^^ff account for collective spin flip processes caused 
by a M0lmer-S0rensen interaction. If this interaction did not exist, the propagator could 
be calculated in the same way as for the driven harmonic oscillator described by ([2]). In 
the limit Vt <^ where riQ{Ji + J3) = 2r]Q'^/6 + O(n^), the time evolution from t=0 to 
t* = 2n/\e\ would create a mapping of quantum states 0(0) (pit*) described by the 
operator 

Ujit*) = expitet*Sl^) 

with 

For m = 2 ions, the operator Uj{t*) performs a conditional phase gate if 6t* = 71/8. For 
C = and weak excitation [Q <^ u), this requires setting the coupling strength Q = Qc 
with 
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In the limit where <^ z/ no longer holds, saturation effects reduce the geometric phase 
$ picked up in the gate operation. For = Qc, we would now have 

()t*^l(l- At) ^ (12) 



SriNtJ 

where Nt = v/\t\ counts the number of trap cycles during the gate operation. For 
?7 = 0.1 and a gate time of 100 trap cycles, $ is reduced by about 7%. The smaller the 
Lamb-Dicke factor gets, the more important saturation effects become for a given gate 
time. The M0lmer-S0rensen interaction contributes a term to the propagator Ui{t*) 
which is now approximately described by 

Ujit*) ^ exp{zet*Sl^) exp{z^t*Sl^) (13) 

with 

For the ratio n/O, 



3^- (15) 

If 7] <C 1, the contribution from the second term oc S"^^ is comparatively small. 

Up to now, we have disregarded the fact that the effective Hamiltonian is valid only 
for times T = ^N, N = 1,2, . . . where 6 = ^{u — e). Therefore, the gate time T needs 
to fulfil |e|T = 2tt as well as 6T = 2ttN , with integer A^. Combining both conditions, 
we find 

e = — - — , N en 
2N +1' 

In writing equation f|T3l) . terms arising from the non- vanishing commutator 
[5*2^^,5'^,^] were neglected. Using the abbreviations VL^ = rjVL{Ji + J3) and ^ms = 
Arj^Vt^ Jq^ / {?)5) , it is convenient to rewrite H^^ff = Ha + Hb with 

2 

since Ha and Hb commute. The time evolution induced by Ha is given by the 
propagator 



UaH) = D(A(t)5,,^) expitmSl^) expit^iSi + Si^)) (16) 
with \{t) = -ze-*2C((]^/e)(e*^* - 1) and ^{t) = (^]„/e)2(et - sin(et)), and for the 



interaction Hamiltonian Hi b = U]^HbUa one finds 

Hi,B = ^{cwms! - si^) + smms., ^.,4)- 

Here, the displacement operator £)(±a;) = C{a)±iS{a) was expressed by the real- valued 
operators C and S. For the special case C = this is equivalent to 

H:,B = ^{D{-AX{t))Sl + D{AX{t))Sl). 



Ion trap quantum gates with amplitude-modulated laser beams 



12 



The last expression shows that the interaction Hamihonian Hi b describes collective 
spin flips between the levels | J.|) and | tt) that go along with displacements of the 
vibrational state. For a phase gate operation, max(|4A(t)|) ~ 2. Minimum uncertainty 
states of motion are not conserved by the interaction. 



3.2. M0lmer-S0rensen gate operation 

The formalism developed so far can be employed to study the M0lmer-S0rensen gate as 
Hamiltonian ([5]) also describes the bichromatic laser field of the M0lmer-S0rensen gate. 
Since the laser frequencies are set close to the blue and red sideband resonance, the only 
difference is that 6 = u — e instead of 6 = ^{u — e), thus changing the values of the 
Fourier components dn used to express D±. Taking into account the leading terms in 
first and second order for the calculation of ([Tj), one finds the effective Hamiltonian 

= -vQiJo + J2)5,,^(ate^(--^) + ae-^(^-C)) _ Jo^Sj, + -^Ji'Sl, (17) 

instead of (ITOl) . Integrating from t = to t* = and neglecting commutators involving 
5*2^^ in the Magnus expansion, yields the propagator 

Ui{t ) ^ exp <^ ( ( ((^0 + J2j + -^Jo 1 Sy^^ ^ Ji S,^^, 

= exp{i\t*Sl^) exp{-i^t* Sl^) (18) 
with 

At* = ((Jo + J2? + ^ (19) 

(20) 

The contribution oc 1/ (26) comes from the counter-rotating term, the red laser coupling 

to the blue sideband and vice versa. For weak excitation, we have 

r. 2r/2fi2 
Ui{t )= exp Y^—^by^^ 

For m = 2 ions, an entangling gate operation is achieved by setting |At*| = | which 
amounts to setting the coupling strength Q to 

= M. (21) 

In the limit where f2 <^ z/ is no longer valid but where \e/ri\ <^ v still holds, we find for 
the correction terms in f|T9l) when keeping f2 = 

|A^1-I^l-7r4^.-%^^V (22) 



4(r/iV,)2 2Nt 



and for the ratio 
u 
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For ?7 = 0.1 and a gate operation that is performed within 100 trap cycles, the correction 
terms to At* have a relative strength of 0.25% and 0.5%, respectively, and the Sl^^ 
interaction is less that 10~^ of the Sy^ term. Therefore, the interaction is quite well 
approximated by using the propagator (fT8|) with /it* set to zero. Then, one obtains for 
arbitrary t 

Ui(t) = b{a{t)Sy^^) exp [i{\t - x sin(et))S'^,^) 

where 

a(t) = ^(Jo + J2)e-*«(e*^*-l) (24) 
e 

A =^(('^o + J2)^ + ^V) (25) 

X =^-^{Jo + J2? (26) 

In the reference frame of the original Hamiltonian ([5]), the laser-ion interaction is 
therefore well described by the propagator 

U{t) = exp{-tF{t)S,)D{a{t)Sy,^) exp {t{Xt - x sin(et))5j,^) . (27) 

This propagator can be used to calculate the dynamics of expectation values of 
interest for the qubits. It is possible to derive simple expressions by tracing over the 
motional states if the vibrational mode is in a thermal state. For this, it is useful 
to note that D{a(t)Sy^^) = J2xbi'^i^)^)P^ where P\ denotes the projector onto the 
subspace spanned by eigenvectors of Sy^^ with eigenvalue A. Moreover, the diagonal 
elements of the displacement operator in the number-state representation are given by 
{n\D{a)\n) = exp(— |ap/2)£n(|ap) where denotes a Laguerre polynomial ^0]. Since 
the generating function of £„(/?) is given by [21] 

n=0 ^ ^ 

summation over a thermal state with number state population p„ = ;^^(^^)" and 
mean phonon number n simply yields 



^Pn{n\D{a)\n) = exp (-\a\^{n + ^) j 



In the case of two ions, C = 0, and an initial qubit state pa = \ \ , the expectation 

value 0(t) = TiQ^Opit)) of the observable O is given by 

Oit) = lTrQ(Ov{(5,' + Si) - 45,e-^l'^l'(^+5) + (5^ - 5^)e-^6l"l'(^+i)}) 
16 

where Trg refers to the trace of the qubit state space and Oy = VOV^ with V{t) = 
exp{—iF{t)Sx) exp{i'y(t)Sy) and jit) = At — xsin(et). As an example, the time evolution 
of (U ii) is explicitely given by 

p^(t) = 1(2 + cos2(2F)) + I cos(2F) cos(47)e-^l"l'("+^) + 1 cos2(2F)e-^^l"l'("+^) 
8 2 8 
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Figure 4. (a) Time evolution of the density matrix elements of two ions in a thermal 
state with n= 2 undergoing a M0lmer-S0rensen interaction with fl = 0.0885 i^,n ~ Q.l, 
e — 0.051/ and C = 0. The calculations are based on ([27]) . The values chosen 
reproduce the curves shown in Figure 3(b) of reference [llj . Counting from above 
at vt = 60, the curves represent the populations Ptt.TT coherences 

/m(pjj^_|-f) and Re{pii^]i). At i>t 250 {t = Air/e), the ions are in a maximally 
entangled state, (b) Same as (a) but with ( ~ n/2. If the gate operation starts in an 
intensity minimum of the amplitude-modulated laser beam, the non-resonant carrier 
oscillations are much stronger. At vt = 250, the quantum state is no longer maximally 
entangled, (c) Fidelity F = {ipmax\p{t)\i^max) of creating the maximally entangled 
state \tp,nax) = (I ii) — i\ TT))/V2 near the optimum calculated from (P7|) . The upper 
curve corresponds to C = 0, the lower one to ^ = 7r/2. The points on top of the upper 
curve represent the fidelity for C = ^^'^ were obtained by a numerical integration 
of the Hamiltonian ([5]) after applying the Lamb-Dicke approximation, (d) Infidelity 
1 — of the gate at vt — 250 for C = and a state with n = 0. The solid line is a 
numerical integration of ([5]) in the Lamb-Dicke approximation, the dash-dotted line 
is based on the full Hamiltonian. The arrow labelled 'a' denotes the optimum Rabi 
frequency predicted by pT|) . '/3' the value of chosen in [11], '7' the Rabi frequency 
predicted by ([25]) . 
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Figure 5. Performance of the (g) gate as a function of the phase C, for a gate 
operation taking place in 25 trap cycles with r] — 0.05, 6 ~ 0.96 and = 0.221 v. 
The figure shows the distance d{Uex, Uid,<p) between the exact propagator Uex obtained 
by numerical integration of (5) and the gate operation Uid,^, = exp{i^Sy^^) (solid line). 
The points on top of the line denote d{Uex, t^(27)) where C/(27) is the propagator given 
by (27), thus demonstrating that this equation is a very good approximation to the 
exact solution. The dashed-dotted line shows d{Uex, Uid,cf,) with Uid,^, = exp(i-|5^) and 
the dotted line d{Uex,Upert) with Upert = exp(i27r?7^fi^/e^5^) as predicted by simple 
second-order perturbation theory. More details regarding the distance measure d are 
given in the text. 

with a{t),'~f{t),F{t) containing the time dependent terms. Other quantities of interest 
could be calculated in the same way. 

A propagator similar to (1271) was calculated in ref. pT] for the case C = 0. The 
authors argued that the non-resonant excitation of the carrier transition could be 
neglected in a first step and obtained in this way a Hamiltonian of the type described by 
([2]) that could be integrated exactly. In a second step, they considered the influence of 
the previously neglected non-resonant excitations. While this treatment yields correct 
results for ^ = 0, it fails to predict the dependence of the gate operation on ( via the 
angle ip as given by fl27|) . Figure H] shows the time evolution of matrix elements for the 
same parameters as used in [TT] for the cases C = ^ind ( = 11/2. In the latter case, 
the amplitude of the non-resonant carrier oscillations is considerable and the input state 
I II) is never perfectly mapped to a maximally entangled state. For C 7^ 0, the effect of 
a non-zero value of ip is fairly small for current gate realizations using hyperfine qubits 
where the Lamb-Dicke parameter r] is considerable. However, it becomes crucial for the 
realization of fast gates on optical qubits with small rj since in this case the gate requires 
a larger value of Q to achieve the same gate speed. 

Figure [5] shows a comparison of the different propagators for the gate operation with 
rj = 0.05 taking place in 25 trap cycles [S = 0.96z/). For the prediction of the required 
coupling strength fl for a gate operation realizing Uid,ii> = exp(z|5'|,^), equation fl25l) 
was iteratively solved to yield At* = n/8. The propagator Uex was obtained from a 
numerical integration of ([5]). Then, Uex was compared to Uid,(f,, to the propagator of 
(1271) . to Uid = exp{i^Sy) and to the prediction Upert = exp(i27rry^f2^/e^S'^) of second- 
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Gate type 




® 


Rabi frequency 




n/u = l/{4:riNt) 


Saturation strength 


7 = 2/(3r7iVi) 


7 = l/(ViV2) 


Coupling ratio 


k/O = 8/(3r/) 





Table 1. Comparison of az (8) Uz and cr^ (g) gates. The first row gives the Rabi 
frequency Q required to peform an entangUng gate operation as a function of rj and 
the gate duration. The latter is expressed as the number of trap osciUation periods 
Nf. The second row hsts the reduction in couphng strength for this kind of gate due 
to saturation effects. The third row compares the unwanted to the desired couphng 
strength. 

order perturbation theory. Since the exact propagator does not perfectly return the 
motional state to the initial state at the end of the gate, the following procedure was 
applied for the calculation of the distance d between the propagators: we assume that 
the ions are initially in the motional ground state and that a cooling mechanisms returns 
the motional state to the ground state at the end of the gate operation without affecting 
the qubit states. This turns the unitary evolution into a quantum process acting only 
on the internal states of the ions. For the comparison of two quantum processes Si, £2, 
the processes are mapped using the Jamiolkowski isomorphism onto density matrices 
pi, p2 for which the distance d{pi,p2) = 1 — i^^{^/^/PlP2^/Pl) is calculated |221[23]. The 
results show that the propagator given by (!27|) correctly predicts the coupling strength 
as well as the operator realized by the gate operation. It also becomes obvious that 
Uid considerably deviates from the operation generated by the Hamiltonian ([5]) unless 

C = o. 

3. 3. Comparison of (Jz ® ciz o-nd (S> o"^ gates 

The main advantage of the cr^ (Jz interaction on optical qubits appears to be its 
insensitivity to changes in the optical path length. In the limit of weak excitation, 
the gate operation tolerates changes that occur within the gate operation as in each 
elementary process a photon is absorbed and another one emitted into the same laser 
beam (this property does not hold for hyperfine qubits since here Raman beams in 
a counter-propagating configuration are used). If higher Rabi frequencies are used, 
the interaction rather becomes az,ip ® Cz^rp which make it susceptible to path-length 
fluctuations within the gate time. Still, if amplitude-shaped pulses are applied (see 
section H]), the gate operation tolerates changes of the path length that occur between 
consecutive applications of the gate. This is not the case for the M0lmer-S0rensen gate 
which becomes robust against changes between gate operations but remains susceptible 
to changes occurring within the gate when the gate is sandwiched between tt/2 pulses 
applied to both qubits. 

The az ® <Jz gate, however, seems to be much less favorable with respect to the 
following criteria: (i) the Rabi frequency that is required for performing the gate 
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operation in a given time, (ii) the strength of saturation effects reducing the couphng for 
the Rabi frequency needed for the gate operation, and (iii) the ratio between the desired 
and the unwanted spin-spin couphngs. Table [1] shows a comparison of the gates with 
respect to these criteria, thus summarizing the results (fTTj) . (fT2|) . (fT5|) for the cr^ ® az 
gate and (12T1) . (|22l) . (l23l) for the a^®a^ gate. For all three criteria, the M0lmer-S0rensen 
performs better. Having a low Rabi frequency is also of interest when it comes to non- 
resonant excitation of other vibrational modes or light shifts induced by excitation of 
far-detuned dipole transitions. 

4. Amplitude-shaped pulses and spin echos 
4.1. Amplitude- shaped laser pulses 

In the limit of fast gate operations, the Hamiltonians ( ITOi) . ( IT71) become sensitive to the 
phase C which is related to the intensity of the bichromatic laser field at the start of the 
gate operation. It is therefore interesting to shape the intensity of the bichromatic 
laser field during the gate operation so that the atomic operator 0{t) = Sj^^(^t), 
with ip(t) = ^^^sm(, appearing in the Hamiltonians becomes time-dependent but 
independent of ( at the beginning and at the end of the gate when the intensity is low. 
In this way, the gate could be made insensitive to ( by an adiabatic process where 0(t) 
evolves from Sj at the start of the gate operations to a (^-dependent operator Sj^^(t) and 
back to 5*^. However, the state a{t) of the vibrational mode generally does not return 
to its original state at the end of the gate under the action of the propagator ([3]) when 
the coupling 7(t) is made time-dependent. There is, however, a class of shaped pulses 
with the property ^(t) =0 that can be constructed in the following way: By applying 
an amplitude-shaped pulse twice with a sign change in the coupling between the two 
pulses, i. e. 72 (^) = — 7i(i), one obtains the propagator 

U = f/_^(2r,r)f/^(r,0) = exp(i2$(r)C2) 

because the first and the second pulse displace the motional state into opposite directions 
but by an equal amount (see ([3])). A quantum state that is displaced along a circle in 
phase space by an off-resonant force of constant magnitude, •yit) = fie*"*, t G [0,27r/e], 
can be viewed as a special case of this pulse form with r = vr/e. For the bichromatic 
gates based on the interactions (|T0|) . f|T7|) . the sign change can be accomplished by 
either shifting the phase ( during the action of the second pulse by an amount 7r/2 (tt), 
respectively, or by changing the overall phase of the laser by vr during the second pulse 

(i.e. n -n). 

Figure [6] illustrates the use of amplitude-shaping in order to make the M0lmer- 
S0rensen gate operation robust against fluctuations in the phase ( between the blue- 
and the red-detuned laser beam. In this example, an entangling gate is accomplished 
within A = 50 trap oscillation periods by a pair of laser pulses with Q(t + r) = —Q{t), 
for t < T with r = vrA/ u. As shown in Figure [6^, the pulses are switched on and off 
within eight trap cycles. After the first pulse, the vibrational state has not returned 
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vf/(2n) Phased 

Figure 6. M0lmer-S0rensen gate operation with two amplitude-shaped laser pulses 
based on the full Hamiltonian ([5]) without Lamb-Dicke approximation. For parameters 
^max = max(|f2|) = 0.167 1', 77 = 0.05, the gate takes place during 50 trap oscillation 
periods. The pulses are switched on and off during 8 trap cycles using a cos^- 
profile. During the second pulse, the phase of the blue- and the red-detuned beam 
is shifted by tt with respect to the first pulse, (a) Time evolution of the populations 
PTT'Pii'PiT +PT1- when starting from state | ]],n = 0) at time t = 0. The dotted line 
shows the coupling strength \il.{t)\/ {2nmax), the dashed line is the average number 
of vibrational quanta, (b) Infidelity of the final state as a function of the phase (. 
The upper curve shows the strong influence of the phase for a gate operation with 
constant coupling strength fl = 0.147 where a high-fidelity operation is achieved 
only for C — 0- For the amplitude-shaped gate, the fidelity is practically independent 
of C- Similar results are also obtained for other input states. For a realistic calculation 
of F, decoherence caused by spontaneous decay of the metastable state would have to 
be taken into account. 

to its initial state. It is only after the second pulse that the correlations between the 
vibrational state and the qubit states vanish again. Using this technique, the initial state 
I t|,n = 0) is mapped to the target state ;^(| It) -|- i\ ii))|n = 0) with an infidelity of 
below 10^^ (see Figure [6)d). This is in sharp contrast to the case of an excitation of the 
same duration with constant amplitude where the infidelity depends on the phase ( and 
varies between 10"'' and 0.2. Similar results are also obtained for other input states. 

4-2. Spin echos 

Spin echo pulses can be combined with amplitude-shaped laser pulses to make the cr^^cr^ 
gate more robust against imperfections. It is possible to implement the conditional phase 
gate operation by having the motional state perform two circles in phase space so that 
the gate pulse can be split up into two separate pulses. Since the quantum states | TT), 
I II) as well as the states | ti), | it) pick up the same phases = = $11, 

it is possible to exchange the states |tt) ^ lii)? lit) ^ iTi) by a collective vr-pulse 
sandwiched between the two gate pulses and to exchange the populations at the end of 
the gate sequence again. The first spin-echo 7r-pulse inverts the direction of the force 
on the motional state so that the motional state returns to the initial state after the 
second spin-dependent pulse. In contrast to the case of shaped pulses without spin-echo, 
there is no need for changing the phase ( of the second pulse or the sign of the coupling 
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strength VL. The spin echo procedure is advantageous for the following reasons: 

(i) The gate becomes more robust against unequal light intensities on the ions. 

(ii) Single qubit phases arising from light shifts are transformed into an unimportant 
global phase. In the context of this gate, light shifts will mainly be due to an 
imbalance in the power of the blue and the red-detuned laser beams and also due 
to very off-resonant excitation of dipole transitions. In addition, a light shift 5is 
occurs if the average frequency 0;^ of the bichromatic light field does not exactly 
coincide with the atomic transition frequency ujq. However, this light shift will be 
fairly small as 5is oc {Vt/SYiujL — ojq). 

(iii) Collective spin flips arising from the term Sy in (fTOl) can be cancelled to first order 
by choosing rotation axes for the 7r-pulses on ion 1 and ion 2 that differ by 90° (x- 
rotation on ion 1 and y- rotation on ion 2). This effectively changes the sign of the 
rotation angle k occurring in ( |T3l) for the second pulse and eliminates the spin flip 
contribution of the interaction. To perform different vr-pulses on both ions requires, 
however, either a different trap frequency that changes the distance between the 
ions by A/4 or an additional laser beam. 

In the limit of short gate operations, spin echos become somewhat less efficient in 
cancelling perturbations described by Sz interactions as the gate interaction oc Sl^.^ 
no longer commutes with for 7^ 0- 

For the M0lmer-S0rensen gate operation, where [Sy^^,Sz] is not a small quantity, 
spin echos seem to be of limited use. If, however, the gate interaction is sandwiched 
between a pair of collective tt/2 pulses to turn it into a (Jz®<Jz interaction, spin echos are 
helpful for cancelling perturbations occurring between consecutive gates. Also, it should 
be noted that a spin-echo like technique was already proposed in ref. jTOj in order to 
cope with number-state dependent ac-Stark shifts that arise if the gate is implemented 
by illuminating ion 1 with a red-detuned laser beam and ion 2 with a blue-detuned laser 
beam instead of using a bichromatic light field for both ions. 

5. Conclusions 

Collective laser-ion interactions with bichromatic laser beams are capable of performing 
both az^Cz gates as well as M0lmer-S0rensen gate operations. The analysis shows that it 
is important to include non-resonant excitation of the carrier transition for the precise 
calculation of the gate operation. While the paper was focused on the case of qubit 
states linked by a weak optical transition, the discussion of the M0lmer-S0rensen gate 
interaction applies also to hyperfine qubits where non-resonant carrier excitation also 
occurs in the limit of fast gate operations. For optical qubits, the required laser beams 
can chosen to be co-propagating which allows for a robust and experimentally easily 
realizable setup where an acousto-optical modulator is used in single-pass configuration 
to create the bichromatic light field. In a direct comparison of the gates, the M0lmer- 
S0rensen interaction seems to be advantageous in terms of required laser power and 
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gate accuracy while the ® (Tz interaction has the advantage of being robust even 
against certain path length fluctuations occuring during the gate operation in the limit 
of weak driving where Sz,^ ~ Sz- For gate durations coming close to T = 2ii/{rji'), 
control of the phase ( between the red- and the blue-detuned laser beams is of vital 
importance unless the gate is performed using amplitude-shaped laser pulses. In this 
case, the requirements are strongly relaxed and the gates appear to be very promising 
for experimental realization. The possibility of using a single laser beam for global 
single qubit and entangling operations opens also interesting perspectives for creating 
multi-particle entangled states with more than two ions. The operations using this beam 
could be combined with an off-resonant strongly focussed beam capable of inducing cr^ 
operations on individual qubits in order to create a larger variety of complex entangled 
states. 
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